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Abstract 

In this paper, we describe direct-sequence code-division multiple-access (DS-CDMA) systems with quadriphase-shift keying 
in which channel estimation, coherent demodulation, and decoding are iteratively performed without the use of any training or 
pilot symbols. An expectation-maximization channel-estimation algorithm for the fading amplitude, phase, and the interference 
power spectral density (PSD) due to the combined interference and thermal noise is proposed for DS-CDMA systems with 
irregular repeat-accumulate codes. After initial estimates of the fading amplitude, phase, and interference PSD are obtained from 
the received symbols, subsequent values of these parameters are iteratively updated by using the soft feedback from the channel 
(3JT)' decoder. The updated estimates are combined with the received symbols and iteratively passed to the decoder. The elimination of 

pilot symbols simplifies the system design and allows either an enhanced information throughput, an improved bit error rate, or 
, greater spectral efficiency. The interference-PSD estimation enables DS-CDMA systems to significantly suppress interference. 

' Index Terms 

Code-division multiple access (CDMA), channel estimation, pilot symbols, expectation-maximization algorithm, iterative 
receiver. 

I. INTRODUCTION 

In mobile communication systems, the wireless channel induces random amplitude and phase variations in the received 
data, with the possible addition of time-varying interference from co-channel users. For this reason, the accuracy of channel 
state information (CSI) at the receiver is critical for coherent detection and demodulation. A number of methods have been 
proposed for estimation of CSI, all of which fall within the broad categories of either pilot-assisted or blind algorithms. Current 
and next-generation cellular protocols such as W-CDMA (Wideband Code Division Multiple Access) and 3GPP LTE (Third 
Generation Partnership Project Long-Term Evolution) specify the use of pilot-assisted channel estimation (PACE) [1|. Pilot 

■ symbols or training sequences are known symbols either multiplexed with or superimposed onto the transmitted data in the time 
\ or frequency domain, with the associated disadvantage of a loss in spectral and/or power efficiency. Moreover, superimposed 

■ PACE is degraded at low signal-to-noise ratios, and multiplexed PACE is unsuitable for fast-fading channels with a coherence 
\ time shorter than the pilot-symbol transmission rate 0, 0. 

t-H . Blind channel- estimation methods offer an alternative approach that avoids the implementation cost of pilot symbols H. 

\ Blind methods typically use second-order statistics of the received symbols for CSI estimation, with shortcomings such as 
• '-j i increased complexity, slow convergence times, and channel-phase ambiguity |5|. In addition, the received interference power 
rS | spectral density (PSD), which is due to both the thermal noise and the time-varying interference, is usually not estimated in 
. the literature spanning both PACE and blind CSI estimation. The accuracy of the interference-PSD estimation is known to have 
a significant impact on turbo-principle (iterative) detection techniques as well as turbo and low-density parity-check (LDPC) 
channel decoding 151, |7l. 

The expectation-maximization (EM) algorithm offers a low-complexity iterative approach to optimal maximum-likelihood 
detection and estimation (§1. l9l. A substantial body of literature can be found on EM-based techniques for data detection, 
multiuser detection, channel estimation, or a combination of the latter. A few representative examples are listed next. A recursive 
estimation of the fading channel amplitude was proposed in [10|. Iterative receivers with EM-based fading-amplitude and data 
estimation using pilot symbols for LDPC-based space-time coding and space-time block-coded orthogonal frequency-division 
multiplexing (OFDM) were studied in [ 1 1 J and [12], respectively. Joint multiuser detection and channel/data estimation for 
uplink code-division multiple access (CDMA) was studied in lfT3l - lfl6l . In flTl . iterative EM estimation and turbo coding 
were studied assuming noncoherent frequency-shift keying modulation and demodulation, which is well-known to be less 
power-efficient than coherent modulation ifTHl . 
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In Ifl9l , an EM estimation approach for turbo-coded single-user iterative CDMA receivers with binary phase-shift keying 
was considered. In ||20l and ll2D . the authors replaced turbo codes with regular LDPC codes; however, Ifl9ll - ||2T1 all featured 
as much as a 9.1% pilot-symbol overhead for channel-amplitude and interference-PSD estimation. Recently, EM-based channel 
and noise estimation techniques were proposed in ll22l and (23 ] for multiple-antenna systems with convolutional coding and 
as much as a 10% pilot-symbol overhead for initial channel estimation. 

Although the primary role of pilot symbols in most cellular standards is channel estimation, pilot symbols often play a 
secondary role in cell, frame, or symbol synchronization. However, alternative methods of synchronization may be used when 
pilot symbols are unavailable [18], l24l . (25). In this paper, a doubly iterative direct-sequence CDMA (DS-CDMA) receiver 
featuring iterative EM channel estimation and iterative detection and decoding without any pilot symbols is presented. The 
general form of the proposed blind channel estimator provides fading-amplitude, phase, and interference-PSD estimates in both 
single-user and multiuser environments, therefore offering an alternative to the methods proposed in ll26ll and ll27ll to rectify 
the phase ambiguity of blind channel estimated The special case of EM channel estimation with perfect phase information 
at the receiver (e.g., by means of a phase-locked loop) is also considered. The proposed iterative receiver is capable of using 
higher-order modulations such as M-PSK and M-ary quadrature amplitude modulation (M-QAM), although quadriphase-shift 
keying (QPSK) is demonstrated in this work for simplicity. In addition, the proposed system uses irregular repeat-accumulate 
(IRA) codes instead of regular LDPC codes for lower complexity Il28l - ll30l . 

The paper is organized as follows. Section II describes the system transmitter and receiver models including coding, 
modulation, and spreading, as well as fading-channel parameters. Section III summarizes the proposed EM-based estimation 
process that uses soft feedback from the channel decoder. Section IV presents the proposed blind method for the initial CSI 
estimation and the possible trade-offs vis-a-vis PACE. Section V shows simulation results, and Section VI offers conclusions. 

A word on notation: lowercase boldface is used to represent vectors, while uppercase boldface represent matrices. E denotes 
the statistical expectation, (-) T is the matrix transpose, * is the complex conjugate, and [a;J is the largest integer smaller than 
x. 
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Fig. 1. DS-CDMA transmitter with QPSK modulation. 

Fig. 1 shows the block diagram of a dual quaternary DS-CDMA transmitter BTI consisting of a channel encoder, QPSK 
modulator, and a direct-sequence spreading generator that multiplies orthogonal chip sequences pr and p/ with the in-phase 
and quadrature modulator inputs. The input to the encoder in Fig. 1 is a binary, independent, identically distributed data block 
of length K, which is denoted by m = [m(l), . . . , m(K)], m(ibit) € [L 0]. 

In 1261 , two different PSK modulations are used on adjacent OFDM subcarriers to resolve the phase ambiguity under slow frequency- selective fading. A 
short pilot sequence is used in [27] to recover the channel phase, making it semi-blind in nature. More importantly, the interference-plus-noise PSD is not 
estimated in \26\ and |27j. 
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A. Encoding, Modulation, and Spreading 

Each 1 x K message vector m is encoded into a 1 x TV codeword b = [6(1), . . . , b(N)] using a systematic, extended IRA 
code J28). IRA codes offer a combination of the linear complexity of turbo encoding and the lower complexity of LDPC 
decoding without compromising on performance. 

The (N,K) IRA code is constructed following the methodology proposed in (291, where the IRA code parameters were 
designed for use on a burst-erasure channel with additive noise, which was shown to be a good surrogate for Rayleigh fading 
channels. IRA codes can be considered to be a subset of low-density parity-check codes and therefore may be represented by 
a Tanner graph [30 1 . Let X(x) — J^." K xl1 an d p{ x ) — Ylt" Pi X% ~ Y represent the variable-node and check-node degree 
distributions of the code's Tanner graph, with (d v ,d c ) being the maximum variable and check node degrees, respectively. Using 
density evolution, for (d v = 8, d c = 7) we obtain the following good choices ||29l : 

A (a:) = 0.00008 + 0.31522x + 0.34085a; 2 + 0.0.06126x 6 
+ 0.28258x 7 

p{x) = 0.62302x 5 + 0.37698x 6 . (1) 

The (TV — K) x TV IRA parity-check matrix can be represented as H = [Hi | H2], where sub-matrix H2 is a ( TV — K) x 
(TV — K) dual-diagonal matrix, and Hi is a randomly-generated (TV — K) x K sparse matrix constructed such that H has the 
degree profile of (Q]). The K x TV systematic generator matrix G is then given by G = [Ir- | H^H^ 7 ] . 

For the simulations in Section V, Gray-labeled QPSK is used with 2 encoded bits mapped into a modulation symbol 
x(k) G {±l,±j} ,k = 1,...,-^. Although QPSK is assumed, the analysis and simulation is easily extended to M-QAM. 
Parallel streams of code bits are each spread using a Gold sequence with spreading factor g chips/code bit before rectangular 
pulse-shaping that produces the real and imaginary components of x(k), i.e., xn(k) = Re (x (k)) and xi (k) = Im (x (k)). In 
practice, an intermediate frequency is used before the carrier frequency upconversion, but the upconversion from baseband to 
the intermediate frequency is omitted for clarity in Fig. 1. 

No channel interleaving is applied to the IRA code due to the inherent interleaving characteristics of the IRA code itself. 
This is because the IRA code can be alternatively represented as a repetition code concatenated with a convolutional encoder 
(accumulator) with an interleaver between them. The interleaver is embedded within the sub-matrix Hi in the Tanner graph 
representation of IRA codes. 



B. Channel Model 

For multiple-access interference (MAI) environments, the channel coefficients are generated using the Jakes correlated fading 
model. The flat-fading assumption is valid when the information bit-rate is low, e.g., 100 kb/s as usually considered in this 
paper, since the multipath delay spread in a typical cellular environment is about 10 /is, which is negligible compared to the 
symbol duration. For completeness, the proposed system and analysis are extended to include frequency-selective channels by 
including multipath components with delays exceeding a chip duration and using Rake receivers fT8l . OTI . as described in 
Section IV-EI Each codeword or frame of TV code bits is divided into two different types of subframes or blocks. One block 
size is set equal to the npB code bits over which the fading amplitude is assumed to be constant. The other block size is set 
equal to nj b code bits over which the interference level is assumed to be constant. 

Each frame comprises TV/2 QPSK code symbols and Ng/2 spreading-sequence chips for each QPSK component. The fading 
coefficient associated with spreading-sequence chip c of either or pj is 

1 Ng n\ 
c = l,...,— (2) 

where E s is the average energy per QPSK symbol, a is the fading amplitude with E [a 2 ] = 1, and tj> is the unknown 
fading-induced channel phase. 



C. Iterative Receiver Structure 

Fig. 2 shows a block diagram of the proposed dual quaternary iterative receiver. The received signal is downconverted, 
passed through chip-matched filters, and despread by a synchronized spreading-sequence generator in each branch, with the 
downconverter and synchronization devices [ 18 1 omitted in Fig. 2 for clarity. Self-interference between the spreading sequences 
of the desired user is negligible because accurate synchronization is assumed at the receiver. Let TVo/2 denote the two-sided 
PSD of the Gaussian noise. For the flat-fading scenario, the complex envelope of the desired user at the k th symbol time with 
active MAI can be written as 

y (k) = C lk/npBi x(k) + n mt (k) + n(k), l<k<- (3) 
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Fig. 2. Iterative DS-CDMA receiver. 



where x(k) is the complex transmitted code symbol of the desired user, n(k) is a complex zero-mean circularly symmetric 
Gaussian noise sample with E |n/c| 2 = iVo, and n mt (k) is the interference at the demodulator due to interfering users lfl8ll . 

El. 

The time-varying MAI is assumed to be generated by interfering users with a structure identical to the desired user, albeit 
the spreading sequences differ and the fading coefficients are independent. The despreading in the receiver tends to whiten the 
interference PSD over the code-symbol passband, and the subsequent filtering tends to produce a residual interference with a 
Gaussian distribution. Thus, the interference PSD due to the combined interference and thermal noise is modeled as additive 
Gaussian noise with a two-sided PSD Iq/2 that is constant over each block of uib code bits but varies from block-to-block. 
This model enables the derivation of an EM estimator for Iq that is used in the demodulator metric and leads to the suppression 
of the interference. 

A receiver iteration is defined as a fixed number of decoder iterations followed by internal EM iterations in the channel 
estimator of Fig. 2, and then a single demodulator metric generation. Let 
i denote the index for the internal EM iteration , i = 1, . . . , i max ; 
j denote the index for the closed-loop receiver iteration, j = 1, . . . ,j max . 



Let 9 



U) 



(i) 



ffi(j) f(J) \ 



represent the estimates of the fading-coefficient and interference-PSD parameters at the i th EM 

iteration during the j th overall receiver iteration. EM iterations commence after the initial channel estimation and decoding, 
which is obtained while the switch in Fig. 2 is set to position 1. The subsequent receiver iterations are performed while the 



switch is set to position 2 in order to refine the initial channel estimate with the aid of soft feedback sjg , j3 
the channel decoder. 



1,2,3,4 from 



III. EM ALGORITHM 

Theoretically, the maximum-likelihood CSI estimator 9 can be obtained from a received data vector y = [y(l), ■ ■ ■ , y(Ni)] 
of Ni code symbols, referred to as the incomplete data, by maximizing the conditional log-likelihood function: 

9 = argmaxln/(y | 9). (4) 



However, the computation of this equation is virtually prohibitive in practice since its complexity increases exponentially with 
the observation window size. In the EM algorithm, the expectation of the conditional log-likelihood of the complete data 
z = (x,y) is iteratively maximized with respect to 9, where expectation is taken with respect to x given y and a previous 
estimate of 9. 

The conditional probability density function (pdf) of z can be written as 

/(z | 9) = /(x, y | 9) = f(y \ x, 0)/(x | 9) = /(y | x, 0)/(x) (5) 
where the last equality is from the independence of the transmitted signal vector x and the CSI parameter 9. Thus, 

ln/(z|0) = ln/(y|x,0)+ln/(x). (6) 
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Since the symbols are independent and circularly symmetric Gaussian noise and interference are assumed, the conditional pdf 

/(y|x,0)is 

ft i m i ( ^ {\y(k)-Cx(k)\*) \ 



fe=i 



Therefore, as |x(A:)| 2 = 1 Vfc, 

1 Nl 

m/(y|x,0) - -iY 1 .l n (J )--^[|y(fc)| 2 + |C'| 2 

-2Re(y*{k)Cx(k))] (8) 

where an irrelevant constant has been dropped. 

E-step: Define the objective function to be the conditional expectation of the conditional log-likelihood of z = (y, x), which 
can be written as 

x{o,0^)=E zly ^\h a f{z\e)] (9) 

where 9^2 is the previous estimate. Using (O and ([8]) and observing that In / (x) in © is independent of 9, and hence 
irrelevant to the maximization, we obtain 



(Mg) = -N/2.hx(I )-jJ^[\y(k)\ 2 +\C\ 2 

k=l 

-2Re (y*(k)Cx^(k))} (10) 



where x^j(k) — ^ [x{k)\ = ^ x | y gu) [x{k)] . Assuming the independence of each transmitted symbol x{k) and the 

independence of x(k) and ($2, and using Bayes' law and the fact that (O can be expressed as a product of N± factors, we 
obtain 

where 



/ (y(k) I x(k) 
f( x (k)\y(k),9§) = ^- 



(*) 



f[y(k)\9 ( ^ 



Pr(x(k)). (12) 



and -|2 



, , , , (13) 

t^o V ■'o 

M-step: Taking the derivative of ( TTOb with respect to the real and imaginary parts of the complex-valued C, and then setting 
the results equal to zero, we obtain the estimate of the fading coefficient at iteration i + 1 as 



1 k=l 

Im (Aw)) = -^E Im (y*w^J(fe)). (is) 



fc=i 

Similarly, maximizing ( fTOb with respect to the interference PSD Iq leads to 

(16) 



k=l 

The fading phase and amplitude can be explicitly estimated from ( TBI and (fT3T >, but that is unnecessary. 

Let Sp , p — 1,2,3,4, be the code-symbol probabilities obtained from the soft outputs of the channel decoder, with 
si = Pr (x(k) = +1) , s 2 = Pr (x(k) = +j) , s 3 = Pr (x(k) = -1) , s 4 = Pr (x(k) = -j). From © and ([121), the expectation 
of x(k) at the i th EM and j th receiver iteration is 

„(i) p0') 1 „-„0') _ p0') _ „•„(•?) n(i) 
0=1 
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where likelihood-ratio i£g r!^ depends on the current CSI estimates as 
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Therefore, for a given receiver iteration, x^(k) and i?^ are updated i m ax number of times using decoder feedback s^. 
the next receiver iteration, after channel re-estimation, the fading-coefficient and interference-PSD estimates are updated, and 
then used at the demodulator and channel decoder to recompute x^ X \k) and Rp~ffi- This process is repeated again for i max 
EM iterations, and the aforementioned cycles continue likewise for subsequent receiver iterations. 

In estimating the fading parameters, we set N\ = in estimating Iq, we choose njB < npB and set N\ = uib/^. 



The EM estimator first finds the value of C^.) for a fading block of size ufb- Then it finds the value of Zw^ for each 



0) 



smaller or equal interference block of size nj b using the value of found for the larger or equal fading block. When pilot 

symbols are used, we set x/)(k) = x{k) for each known pilot bit, and there are no EM iterations if only known pilot bits 
are processed in calculating the channel estimates. The application of the EM algorithm is to obtain both channel-coefficient 
and interference-PSD estimates, which differs from lfTT1 - lfT3l where the emphasis is on data detection, and noise statistics are 
assumed to be perfectly known. 

Let I = 1,2 denote the two bits of a QPSK symbol, and t>i,va denote the corresponding log-likelihood ratios that are fed 
back by the channel decoder. From l2~Tl and |[32~1 Eqn. 6], the demodulation metrics (extrinsic information) zf (k),l = 1,2 
for bits 1, 2 of symbol k that are applied to the channel decoder are shown at the top of the next page. 

The number of EM iterations and the receiver latency are reduced by applying a stopping criterion. Iterations stop once 
C$ is within a specified fraction of its value at the end of the previous iteration or a specified maximum number is reached. 
The fraction should be sufficiently small (perhaps 10%) that the performance loss will be insignificant. 



o,(i) 



IV. BLIND CSI ESTIMATION 

The EM algorithm in Section III generates updated CSI estimates as shown in (TT4T> — <TT~6b after the initial coherent demodulation 
and decoding of receiver iteration j = 0. In [19|-[21|, the initial CSI estimates were obtained with the aid of pilot symbols. 
In this section, two methods for blind estimation of the initial CSI parameters 0^ j = )>^o(] )J are presented, 

with the special case of perfect phase information at the receiver examined first. 



A. Perfect Phase Information at Receiver 

The carrier synchronization provided by a phase-locked loop in several second and third-generation cellular standards such as 
IS-95 and CDMA20001x can be exploited to obviate the need to estimate the channel phase (which is also potentially provided 
by 2% piloting [27]). Assuming perfect phase information at the receiver, the fading amplitude is real-valued and nonnegative, 
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and (fT3T > does not have to be computed. A simple heuristic estimate (denoted as blind method I) of ) , 

each fading block can be obtained from the received symbols as 

ufb/2 



f(0) 
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for 



(21) 



(22) 



(23) 



, j is the difference between that power and the 



represents the average power of the received symbols, and D — ^ 

estimated average power of a desired symbol. Equation ( |2"TT i would provide a perfect estimate in the absence of noise and 
interference. The parameter ft, > is chosen such that (^O^ j j does not exceed some maximum value. Ideally, 
h is a function of E s /Nq, but in this paper a constant h = 0.1 is always used for simplicity. 



B. Complexity Analysis 

Although the EM estimation is a relatively low-complexity iterative approach to maximum-likelihood estimation, it consumes 
a much larger number of floating-point operations than pilot-assisted schemes do. To evaluate the complexity of the EM estimator 
in terms of required real additions and multiplications per block of N\ code symbols, each complex addition is equated to two 
real additions, each complex multiplication is equated to four real multiplications, and divisions are equated with multiplications. 
Equations (fl4]>— (fT6b require j m ax*max (6iVi + 4) real additions and j max i m ax. (12iVi + 4) real multiplications. Equations (18) 
and (19) require 6j max imax real additions, 30j max «max real multiplications, and the computation of 4 exponentials. The initial 
estimates calculated using (1211— (1231, which only need to be computed once prior to the first EM iterations, require 2N\ real 
additions, 8A r i + 7 real multiplications, and the computation of the maximum of two real numbers. A PACE receiver that uses 
only pilot symbols for CSI estimation requires 6Ni + 4 real multiplications and 12A^ + 4 real multiplications to compute 
(TBI— (fToll once and does not need to compute the other equations. Thus, EM estimation increases the amount of computation 
for CSI estimation by a factor of more than j max i max relative to PACE. 



C. No Phase Information at Receiver 

The initial CSI estimates proposed in (|2TT > and (l22l for blind method I are expected to be degraded significantly when the 
phase information is also unknown, since an arbitrary initial phase value (e.g., radians) must be assumed. To circumvent this 
problem, the initial receiver iteration consists of hard-decision demodulation and channel decoding, after which each decoded 
bit is used as x^P) Jk) in (fT~4T> — (fT~6l>. This step is followed by the regular EM estimation process in subsequent receiver 
iterations. This approach for the initial CSI estimates, which is referred to as blind method II in the sequel, results in increased 
receiver latency relative to the previous method when phase information is not available. 



D. Blind-PACE Estimation Tradeoffs 

The previously proposed iterative DS-CDMA receiver with PACE 1 19|-[21 1 is considered as the benchmark for comparison 
with the proposed receiver. Assuming an identical transmit-power constraint and information bit-rate in both cases, the 
elimination of pilots creates the following possibilities for methods I and II: 

• (Case A) An increase in the number of transmitted information symbols. 

• (Case B) An increase in transmitted information-symbol duration. 

• (Case C) An increase in the number of transmitted parity symbols (lowered IRA code rate). 

The modifications listed above offset the loss in system performance due to the degraded CSI estimation obtained from 
blind methods I and II with respect to PACE. The no-pilot cases A, B, and C have the same transmitted frame duration as 
the frame with pilot symbols. Cases A, B, and C provide the most favorable throughput, spectral efficiency, and bit error 
rate, respectively. Numerical evaluations of each of these cases are presented in the next section. Although a correlated fading 
model is assumed in the simulations, no filtering is used to exploit this correlation in order to maintain the robustness of the 
proposed estimator. 
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V. SIMULATION RESULTS 

In all the simulations, the block sizes are equal, and the information-bit rate is 100 kb/s. Increasing the block sizes increases 
the accuracy of the EM estimators, but decreasing the block sizes allows closer tracking of the channel parameters and 
includes more diversity in the receiver computations. In most of the simulations, except where stated, we set riFi=nFB = 40 
and spreading factor g = 31. The number of closed-loop receiver iterations is set to j ma x = 9, as there is insignificant 
performance improvement for j max > 9. The number of internal EM iterations is i max = 10. There is one decoder iteration 
per receiver iteration. A IRA code (data block size K = 1000) with sum-product algorithm decoding [18] is used without 
channel interleaving. The IRA code is rate-1/2 when PACE is used. Jakes correlated fading of the desired signal and a mobile 
velocity of 120 km/hr are assumed. Flat fading is assumed in most of the simulations, whereas a frequency-selective channel 
is examined in Section fV-EI The iterative PACE receiver considered for comparison contains 9.1% pilot-symbol overhead, 
which has been shown to have a decoding performance close to the conventional 3GPP LTE receiver lED . For each scenario 
tested, 5000 Monte Carlo simulation trials were conducted. To avoid repetition, a selection of representative examples out 
of the many possible combinations of channel coding, phase information, interference models, and no-pilot modifications are 
presented next. 

The bit error rate (BER) is calculated as a function of E^/No, where Eb = (N/2K)E S is the energy per bit. The information 
throughput is a vital performance criterion in addition to the BER. One of the primary motivations in removing pilot symbols 
is the expectation of achieving greater throughput, even though the BER performance may be degraded marginally. We define 
throughput R as 

information bits in a codeword 

R = X (1 - BER) bits/s. (24) 

codeword duration 

A. Single-user environment, perfect phase knowledge 




Fig. 3. BER versus Et/No for IRA-coded iterative receiver in single-user environment with phase provided by PLL. 
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Fig. 4. Information throughput versus E^/Nq for IRA-coded iterative receiver in single-user environment with phase provided by PLL. 



For the first set of results in Figs. 0431 a single-user environment and perfect phase knowledge at the receiver are assumed. 
Fig. [3] displays the BER versus Eb/Ng for an IRA-coded iterative receiver operating with perfect CSI, PACE, blind method I 
with cases A, B, and C, and blind method II with cases A and C, respectively. The key observation is that blind method II 
is worse than method I by 2 dB at BER = 1CP 3 for both case A and case C, which illustrates the well-known sensitivity of 
the EM algorithm to the accuracy of the initial estimates. 

The addition of extra parity bits to blind method I (case C, rate- 1000/2200) offers the greatest improvement in BER, 
surpassing even the rate-1/2 code with perfect CSI at high E^/Nq. The increase in number of information symbols (case 
A) results in the worst BER performance with a separation of 1 dB and 0.5 dB from PACE and case B at BER = 10~ 3 , 
respectively. 

The various scenarios featured in Fig. [3] were also tested under a slow-fading channel with mobile velocity of 10 km/hr. It 
was observed that all the BER curves were shifted towards the right by up to 7 dB at BER = 10~ 3 , but the overall trends 
among the different cases remained the same. 

Fig. |4]exhibits information throughput R versus Eb/N n for the IRA-coded iterative receiver with the scenarios of Fig. [3] The 
throughput advantage of case A is achieved even though no pilot symbols are used at all; i.e., the initial estimation is totally 
blind. It is evident that increasing the symbol duration or adding additional parity information does not give the proposed blind 
methods any significant advantage in throughput over PACE. Both blind methods with cases B, C and PACE provide about 
20% less throughput than the receiver with perfect CSI. 



B. Multiuser environment, unknown phase 

A 4-user interference environment with equal mean bit energies for all users at the receiver, E^/Nq — 20 dB, and no 
phase information at the receiver is examined next. It is assumed that both the interference levels and the unknown phase are 
constant during each subframe. Each interference signal experiences independent Jakes correlated fading and uses independent 
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data and Gold sequences with respect to the desired signal. The simulation uses chip-synchronous interference signals, which 
is a pessimistic worst-case assumption (3T). Two variations of CSI estimation are examined here: partially adaptive with 
only fading coefficient being estimated using (fl4] >. ([T5l l. and I^X set equal to No for all subframes; and fully adaptive 

estimation of both c9.) and PA\ using (flit . ( fl3T >. and ( fT6] l. 
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Fig. 5. BER versus E^/Nq for IRA-coded iterative receiver affected by MAI from 4 users, fully and partially adaptive estimation, and unknown phase. 

Fig. [5] displays IRA-coded BER versus Eb/No for partially and fully adaptive CSI estimation per fading block and case C 
for both blind methods. The mismatch of Io and the true value of Iq at the demodulator and decoder results in a high error 
floor for the partially adaptive cases. The intuition behind the error floor is that the partially adaptive estimator overestimates 
the true signal-to-interference-plus-noise ratio (SINR) by disregarding the MAI, with the degree of overestimation increasing 
with SINR. For IRA codes, it was shown in [33] that both under- and overestimation of the SINR degrades the IRA decoder 
performance. The fully adaptive estimation offers a more accurate SINR estimate and, hence, suppresses interference and 
reduces the error floor significantly. This interference suppression is achieved without using the far more elaborate multiuser 
and signal cancellation methods that could be implemented in a DS-CDMA receiver. For both partially and fully adaptive 
estimation, it is observed that blind method II now outperforms method I due to better phase estimation, whereas both blind 
methods outperform PACE at BER = 10~ 3 due to the added parity information. 

Fig. [6] demonstrates the IRA-coded receiver throughput offered by the proposed methods under MAI from 4 users. The blind 
methods always provide a better throughput compared with PACE; for example, method I with case A is superior by 9% to 
both PACE scenarios when E^/Nq > 5 dB. It is observed that both partial and fully-adaptive estimation methods offer a similar 
asymptotic throughput, which indicates that partial CSI estimation may be sufficient for applications with a non-stringent BER 
criterion. On the other hand, error-critical applications requiring less than BER = 10~ 3 must use the fully adaptive CSI 
estimation, as seen from Fig. [5] 
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Fig. 6. Information throughput versus E^/No for IRA-coded iterative receiver affected by MAI from 4 users, fully and partially adaptive estimation, and 
unknown phase. 



C. Varying fading-block size, unknown phase 

In urban mobile environments, the phase can be expected to change significantly after approximately s to 2^4 s, where 

J d Id 

fd is the maximum Doppler shift. For the assumed mobile velocity of 120 km/hr, this time range corresponds to roughly 10 to 
40 code bits at 100 kb/s. The fading and interference block sizes nps = niB are therefore varied accordingly, and no phase 
information is assumed to be available at the receiver for the next set of results. 

Fig. [7] displays fully adaptive IRA-coded BER versus E^/Nq for blind methods I and II with case C, 9.1 % PACE, and 
perfect CSI decoding for ufb = 10 and 40 in a single-user environment. An improvement of 1 to 2 dB was observed for all 
methods for the smaller fading-block size of nps = 10 due to the increased fading diversity. The throughput with case A is 
shown in Fig. [8] It is observed that the throughput gains of the proposed blind methods over PACE (roughly 9% at intermediate 
to high Eb/No) are preserved even when the phase is initially unknown at the receiver. 

D. Varying MAI, unknown phase 

IRA-coded iterative receiver performance with blind method II, case C is examined for 3 and 6 MAI signals with equal mean 
bit energies for all users at the receiver in Fig. [9] The partially adaptive estimation is unable to cope with the interference caused 
by 6 MAI signals regardless of the spreading factor, whereas the fully adaptive estimation offers a substantial improvement in 
BER. The benefit of an increased spreading factor (g — 127 versus g = 31) is more apparent at low bit error rates for fully 
adaptive estimation. For example, the fully adaptive estimation with 3 MAI signals improves by a factor of approximately 5 
dB at BER = 10 -5 , despite nonorthogonal spreading sequences and imperfect CSI. 
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Fig. 7. BER versus E^/Nq for IRA-coded iterative receiver in single-user environment, varying rips, and unknown phase. 



E. Multipath channel 

A DS-CDMA system can exploit a frequency-selective fading channel by using a Rake receiver. As an example, we assume a 
channel with three resolvable multipath components (with known delays) of the desired signal and a Rake combiner with three 
corresponding fingers. The multipath components undergo independent fading across the fingers, but follow the Jakes correlated 
fading assumption over time. The multipath components follow an exponentially decaying power profile across the fingers, 
i.e., E [ai] 2 = e _ ' i_1 \ I = 1, 2, 3. Each interference signal has the same power level in each finger and undergoes independent 
Jakes correlated fading. The assumption of independent multipath fading amplitude and phase coefficients for the desired signal 
allows us to apply the proposed EM-based channel estimation scheme separately in each finger. The Rake combiner performs 
maximal-ratio combining (MRC) of the received symbol copies based on channel and interference-PSD estimates computed 
at all fingers. The MRC decision statistic obtained from the Rake combiner is then passed to the QPSK demodulator metric 
generator, which generates soft inputs for the common channel decoder. The channel decoder soft outputs are fed back to the 
three channel estimator blocks, which then recompute updated channel coefficients, as described in Section III. 

Fig. [10] displays the Rake receiver performance for various levels of MAI with Method II under case C, where all users 
have length- 127 Gold sequences. It is observed that the additional diversity due to Rake combining improves performance as 
expected, but the performance disparity between partially and fully adaptive estimation remains large. 

VI. Conclusions 

It has been shown that pilot symbols are not essential to the effectiveness of DS-CDMA receivers with coding, coherent 
detection, and channel estimation. If the pilot symbols are replaced by information symbols, the throughput increases relative 
to PACE whether or not interference is present. If the BER is the primary performance criterion, then replacing the pilot 
symbols by parity symbols gives a lower BER than PACE. If the spectral efficiency is of primary importance, then extending 
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Fig. 8. Information thi'oughput versus E^/Nq for IRA-coded iterative receiver in single-user environment, varying npB, and unknown phase. 



the symbol duration after the removal of the pilot symbols offers an improvement relative to PACE, albeit at the cost of a 
slight increase in the BER. 

The estimation of the interference PSD has been shown to enable the significant suppression of interference. This suppression 
is achieved without using the far more elaborate multiuser and signal cancellation methods that could be implemented in a 
DS-CDMA receiver. 
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